Quantum phase transitions in three-leg spin tubes 
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We investigate the properties of a three-leg quantum spin tube using several techniques such as 
the density matrix renormalization group method (DMRG), strong coupling approaches and the non 
linear sigma model (NLctM). For integer S, the model proves to exhibit a particularly rich phase 
diagram consisting of an ensemble of 2S phase transitions. They can be accurately identified by the 
behavior of a non local string order parameter associated to the breaking of a hidden symmetry in the 
Hamiltonian. The nature of these transitions are further elucidated within the different approaches. 
We carry a detailed DMRG analysis in the specific cases S — 1. The numerical data confirm the 
existence of two Haldane phases with broken hidden symmetry separated by a trivial singlet state. 
The study of the gap and of the von Neumann entropy suggest a first order phase transition but at 
the close proximity of a tricritical point separating a gapless and a first order transition line in the 
phase diagram of the quantum spin tube. 

PACS numbers: 75.10.Jm, 75.10.Pq, 64.70. Tg 
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I. INTRODUCTION 

Frustrated spin models in one dimension have attrac- 
ted attention for both the uniqueness of their characte- 
ristics and the diversity of their properties. In contrast 
to higher dimensional spin systems, quantum spin chains 
have no long range order. If there is no frustration, the 
properties of the chain are essentially governed by the 
parity of the spin : the Heisenberg spin chain for ins- 
tance has a gapless spectrum and algebraic correlations 
when the value of the spin is a half-integer whereas it has 
a gap and exponentially decaying correlations when the 
spin is an integer 1 . When frustration is present, the pro- 
blem gets much more complex and the possibilities for the 
low-energy spectrum are also broadened. An illustrative 
example is given by the spin ladder with S — 1/2 and ad- 
ditional diagonal couplings. Depending on the strength of 
the frustrating couplings, the ground state of the system 
can be described in terms of rung singlets, short-range 
valence bonds 2 , or would eventually dimerize 3 -'. The tran- 
sitions between some of these phases have been proposed 
to be deconfmcd quantum critical points which could sup- 
port fractionalized spinona^. 

Another family of problems concerns the integer spin 
ladders. The comprehension of integer spin chains have 
considerably improved since the discovery of the AKLT 
Hamiltonians 5 and the early work of den Nijs and Rom- 
melse^. In particular, the ground state of the spin-1 Hei- 
senberg chain is now well understood : it displays a subtle 
hidden topological degeneracy-^, associated to a non- 
vanishing non-local string order parameter 6 and supports 
edge states. The question of the preservation of the topo- 
logical order when couplings between different chains are 
introduced is an open issue. It is believed that this order 
should be highly sensitive to perturbations. As a mat- 



ter of fact, a simple coupling between two spin-1 chains 
leads rapidly to the destruction of the topological ordeiM^, 
reflecting the fragility of the edge states towards the per- 
turbation (see also Ref.llll) However, it is also possible 
to maintain the topological phase by adding frustrating 
nearest-neighbor interactions. In this case, a direct first- 
order transition between two different topological phases 
can be observed^. The question of the stability of the 
topological order in spin ladders is of crucial importance 
if one thinks of these systems as intermediates between 
Id and 2d systems and regards them as a pathway to dis- 
cover a spin liquid behavior in two-dimensional systems. 

In this work, we investigate the presence and nature of 
topological phases in an asymmetric three-leg quantum 
spin tube with integer spin quantum numbers. The tri- 
angular spin tube has already been extensively studied 
in the spin-1/2 case. Abelian bosonisation techniques^ 
arguments suggest that the system is gapped when the 
tube is symmetric and maximally frustrated. It is interes- 
ting to introduce asymmetry among the coupling in each 
triangle. The model with the asymmetry has been stu- 
died by density-matrix-renormalization-group (DMRG) 
algorithm. Recent DMRG calculations^ have demons- 
trated that the dimer order is unstable against a small 
but non zero anisotropy coupling, that eventually drives 
the system into a critical phase. 

Much less is known in the case of integer spins. The 
triangular geometry provides a simple and natural way 
to introduce frustration, and we thus hope to find uncon- 
ventional behaviors. Here, one coupling between two legs 
is varied, thus controlling the strength of the frustration 
(see Figure [l]) in order to explore a large phase diagram. 
The possibility of quantum phase transitions with decon- 
fined spinons is also an interesting question. 

Besides DMRG, a group of methods that can be used 



to investigate this problem are the large-5 approaches. 
Among them is the non linear sigma model (NLctM) 
which furnishes crucial informati on regarding the spec- 
trum of spin chains and ladders^^. Spin models with 
triangular geometry are described in the continuum by 
a 50(3) rotation matrix field, in contrast to collincar 
antiferromagnets for which the NLctM theory involves a 
single unit vector field 16 . 50(3) NLctM are characterized 
by the absence of a topological term in the actional an(1; 
in d > 1, by a non-trivial fixed point with an enlarged 
50(4) symmetry^]. Even without topological term, in- 
teger and half-integer spin beha ve diff erently due to the 
occurrence of topological defect ^^21. It remains to be 
seen how this scheme is perturbed by the introduction of 
an anisotropy in the triangular geometry. 

We determine the phase diagram of the anisotropic 
spin tube with integer spin 5 by gathering together the 
results obtained from diverse methods : strong coupling 
expansion, large-5 approaches and DMRG. We find that 
the tube supports 25 quantum phase transitions when 
the anisotropic coupling is varied. The nature of the tran- 
sitions is debated. We begin in section [TT] with the proper 
definition of the model and introduce its strong coupling 
limit. Different phases are delimitated depending on the 
value of the quantum spin J of each triangle. In section 
|III[ we develop the notion of string order parameter and 
we show how the spin tube model can be rewritten in 
terms of a local Hamiltonian with a discrete Z 2 x Z2 
symmetry. This hidden symmetry is broken when J is 
odd and remains unbroken when J is even. To unders- 
tand the nature of the phase transition, we turn in the 
third part to the large-5 approaches. We begin with a 
spin-wave analysis to determine the low-energy modes of 
the model. We then derive the NLuM and the associated 
Renormalization Group (RG) equations. In our deriva- 
tion, we put a careful emphasis on the evaluation of the 
total Berry phase of the tube. We find 25 special values 
of the anisotropic coupling corresponding to a non tri- 
vial Berry phase. Then, we focus on the special case of 
the spin-1 tube with a strong coupling approach and a 
DMRG study. The DMRG results reveal the presence of 
two quantum phase transition points, in adequacy with 
the predictions of the strong coupling limit. The order of 
the transition is proposed to be first order but the nu- 
merical data also strongly suggest the proximity of the 
system to a tricritical point. Finally, we provide a nu- 
merical phase diagram for the spin-2 tube where various 
even/odd J phases compete. 




Figure 1: (color online) The three different parameters defi- 
ning the coupling of the three-leg spin tube. 



nian reads : 

H = Fy + H ± (1) 

^11 = "^11 Z_^i (^. a ' Sj+l,a) 

H± = J±_ 2_^ (Si,3 • Sj ; i -I- Si : 2 • Sj,3 4- aSi t i • Si^) > 

i 

with i = 1,...,N being the intra-chain index and a — 
1, 2, 3 being the rung index. The point a = corresponds 
to the unfrustrated open ladder while a = 1 is also spe- 
cial because of its translation symmetry in the transverse 
direction. 

B. The classical case 

We start by determining the classical configurations of 
spins which minimize the energy of each triangle by re- 
placing the spin operators S a with classical vectors 5n a . 

For a > 0.5 the solutions that minimize the energy are 
of the kind of the coplanar solution of Figure [2] (b) . 
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II. THE MODEL AND SOME SIMPLE LIMITS 

A. The model 

The anisotropic triangular spin tube is a quantum lad- 
der problem defined by three relevant parameters (Figure 
fll) : the parallel coupling J» , the perpendicular coupling 
J± and the anisotropy parameter < a. The Hamilto- 
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In the extreme limit a — > 00, the two vectors ri! and 
ri2 point in opposite direction and the third spin is es- 
sentially free. The system reduces then to the problem of 
one single chain. On the opposite, decreasing a one enters 
the regime < a < 0.5 in which the lowest energy state 
is an alternated collinear configuration of Figure [2ja). 




states | J, 5*12) between a = 1 and a = is 



Figure 2: (color online) Top : The collinear configuration 
which minimizes the energy for < a < 0.5. Bottom : the 
coplanar configuration which minimizes the energy for 0.5 < 

Q. 
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The first level crossing happens for a = jr-g- The last 
level crossing occurs at a = 0.5. From this last result we 
can conclude that both, classically and quantum mecha- 
nically, the point a = 0.5 corresponds to the entrance 
into the unfrustrated open ladder regime given by a = 0. 
On the other side of the isotropic point, a > 1, the se- 
quence of ground state is given by : 



|0, 5) -*• |1, S - 1) -> ... -> \S- 1, 1) -> |5, 0) (5) 



The first crossing takes place at a = ^M- and the last 
one occurs at a = 1 + 5. After this point, the triangle 
consists of two spins coupled into a singlet and an isolated 
spin. For instance, for 5=1, there is a level crossing at 
a = 0.5 between the singlet state |0, 1) and the triplet 
1 1,2) and another one at a = 2 between the singlet and 
the triplet 1 1 , 0) . On Figure |3j we plot the evolution of 
the main energy levels for 5 = 2. 



In this regime the physics becomes the one of an open 
unfrustrated ladder. 

Note that the collinear state and |2| are both conti- 
nuously degenerate but have a different degree of dege- 
neracy. For < a < 0.5, any alternated collinear confi- 
guration minimizes the energy. Thus, choosing a ground 
state is equivalent to picking up an oriented axis. For 
0.5 < a, all the classical ground states are given by a 
global rotation of the triad (n 1 ,n 2 ,n 3 ). This, in turn, 
requires to choose an oriented axis and an angle. 




C. Quantum spins : the decoupled limit 

Introducing the triangle spin J = Si 4- S2 + S3 and the 
bond spin S 12 = S x + S 2 , the rung Hamiltonian reads : 



Figure 3: (color online) Spectrum of a single triangle for 
S = 2 as a function of a as given by the two sequences Q 
and (TSu . The higher energy levels are not represented. 
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where we have replaced the spin operators by their eigen- 
values. To determine the ground state, we need to label 
each state by their value of total spin J and their inter- 
mediate spin 5 12 . For a = 1, the 5i 2 levels are degenerate 
and the ground state is obtained for the smallest value 
of J. Thus, the ground state is the singlet state |0, 5) 
(if J = 0, 5i 2 = 5 necessarily). When turning on the 
anisotropy, other levels will compete with this state. It 
is straightforward to show that the sequence of ground 



Thus, in the strong rung coupling limit Jn = 0, there 
are 25 transition points in < a < CO. If we now add 
a small longitudinal coupling J11 <ti J±, we can still ex- 
pect that 25+1 different phases are present. However, 
we need to know how to effectively make the distinction 
between them. As we will see in the next chapter, this 
can be achieved with a non local string order parameter. 
The nature of the phases corresponding to J — and 
J^0 are clearly different. In the former case, the tube 
just consists of a trivial superposition of singlets. We will 
refer to this phase as the singlet phase. In the latter, 
the properties of the tube are more similar to those of a 
single chain with 5 = J spins, and we will refer to them 
as Haldane-like phases. 



III. INTEGER S CASE AND HIDDEN 

SYMMETRY 

A. Hidden symmetry and string order parameters 

In order to characterize the different phases suggested 
by the strong coupling analysis, we would like to find 
a suitable order parameter enabling us to describe the 
phase transitions. Usually, different phases are characte- 
rized by (local) order parameters, which detect sponta- 
neous symmetry breakings. However, in some cases this 
standard approach does not work. This includes, in par- 
ticular, the Haldane phase of S = 1 chain : it has no 
local order parameter but still is a distinct phase separa- 
ted from a trivial phase by a quantum phase transition. 
In order to characterize the Haldane phase, the non-local 
"string order parameters'^, one of which is 
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(6) 



is useful. It has been confirmed that it is non-vanishing 
within the Haldane phase but is zero in a trivial phase 
(for example the large anisotropy phase of Ref. [S]) . 

The problem with the non-local order parameter such 
as Eq. ^ , in general, is that it is not quite clear if there is 
necessarily a phase transition between two states, when 
a non-local order parameter vanishes in one state but is 
non-zero in the other. Kennedy and Tasaki^ clarified the 
meaning of the string order parameter Mn, as an order 
parameter measuring a spontaneous breaking of hidden 
discrete symmetry. Namely, there exists a non-local uni- 
tary transformation which transforms the Hamiltonian 
% to a Hamiltonian % with short-range interaction and 
with a discrete Z2 x Z2 symmetry. The Z2 x Z2 symmetry 
is hidden in a non-local way in the original Hamiltonian 
H. 

The string order parameter (fGJ) is transformed by the 
same non-local unitary transformation to the standard 
ferromagnetic order parameter. Thus, non- vanishing 
string order parameter ([6]) implies a spontaneous brea- 
king of the hidden Z2 x Z2 symmetry. The spontaneous 
symmetry breaking clearly distinguishes the phases. In 
this sense, the string order parameter indeed qualifies as 
an order parameter, despite its nonlocality. 

The hidden Z2 x Z2 symmetry breaking also implies 
4-fold groundstate degeneracy. This appears contradic- 
tory to the uniqueness of the groundstate in the Hal- 
dane phase. However, the non-local unitary transforma- 
tion only works for the open boundary condition. Thus 
the hidden Z2 x Z2 symmetry breaking implies 4-fold 
groundstate degeneracy of the original Hamiltonian 1-L 
only in the open boundary condition. This degeneracy ac- 
tually corresponds to the existence of the edge states with 
spin Sb = 1/2 at both ends. 

The appearance of the edge states can be unders- 
toocP"^ in the Valence Bond Solid (VBS) picture where 
each spin 1 is seen as a triplet of spins 1/2. Spin 1 at each 



site is first decomposed into two spin 1/2's. Each consti- 
tuent spin 1/2 is then coupled to a spin 1/2 in the neigh- 
boring site to form a singlet. This would give a simple 
dimerized state of a spin 1/2 chain. However, projection 
to the triplet sector within each site gives a nontrival 
state for S — 1 chain. If we consider this state on a fi- 
nite chain with the open boundary condition, unpaired 
spin 1/2 is left free at each end. Namely, spin 1/2 degree 
of freedom appears at the ends. The above construction 
actually gives exact groundstates for a special, solvable 
Hamiltonian. For other models, the constructed state is 
of course not an exact groundstate. However, the ap- 
pearance of the edge states is a common feature within 
the S — 1 Haldane phase. As a consequence of the edge 
states, the groundstates of an open chain is asymptoti- 
cally 4-fold degenerate. As mentioned above, this corres- 
ponds to the spontaneous breaking of the hidden Z2 x Z2 
symmetry. 

Thus, the hidden Z2 x Z2 symmetry breaking characte- 
rizes the S ~ 1 Haldane phase, unifying the string order 
parameter and the edge states. However, it should be also 
noted that this picture is only valid in the presence of the 
global Z2 x Z2 symmetry. We will, in Sec. |VIH| discuss 
from the perspective of rec ent, m ore general characteri- 
zation of the Haldane phas e 1 21 1 22 - 1 . 

Now let us move on to our problem of the spin tube 
with integer spin. Naturally, ladders/tubes are more com- 
plicated than the single chain, and various generalizations 
of the string order parameter have been proposed. Howe- 
ver, as we have discussed for the single chain, generally 
there is no guarantee that a non-local "order parameter" 
really qualifies as an order parameter. Therefore, in this 
paper, we first generalize the hidden Z 2 x Z 2 symmetry 
to the tube. Then we identify the corresponding string 
order parameters, which detect spontaneous breaking of 
the hidden Z2 x Z2 symmetry. 

Following Ref. [231 the Kennedy- Tasaki transformation 
generalized to the tube can be written as : 



V = [] exp {i*3*3l) 



(7) 



with Jj = S;.i + Si,2 + Si,3. We impose the open boundary 
condition on the tube (along the leg direction). 

It is straightforward to show that the spin operators 
transform into : 
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The natural generalization that comes to mind is to 



define the two string order parameters : 

(^) = |Jf n oo ^fe X p^t/r)^) (8) 

(0-> = | ^ oo ^«p(«rgj?)^ (9) 

Applying the unitary transformation ufh, they reduce to 
the local ferromagnetic order parameters : 

(6 a ) = {vo a v- 1 ) = lim (j?j£), (io) 

for a = x, z. Now, let us consider the Hamiltonian. This 
transforms into : 



H = H 
H 
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Note that the rung part is invariant under the non-local 
transformation. The new Hamiltonian still consists of lo- 
cal interactions but the global continuous SU(2) symme- 
try of the original Hamiltonian has been hidden and only 
a discrete Z 2 x Z 2 symmetry remains explicit : it is now 
only invariant under the rotation of all spins around the 
x and 2 axis by an angle of 7r. 

We suggest that this "hidden" (non-local) symme- 
try and its associate string order parameters delineate 
the phases of the system. In the strong coupling limit 
(J|| <C Jj_), the phase diagram of the spin tube consists 
of 25 + 1 phases, labelled by the spin index J, analogous 
to the Haldane state for a spin- J chain. It has been de- 
monstrated by one of us^ that not all Heisenberg spin 
chains, but only the ones with J odd, do break the hidden 
Z 2 x Z 2 symmetry and possess a non zero string order 
parameter. Thus, as the anisotropy parameter a is va- 
ried in the spin tube (with Jj_ 3> Jn), we will encounter 
a succession of phases with the string-order parameters 
(pH)- ([9]) alternatively vanishing and non vanishing. 

It is also interesting to consider a disorder parameter 
which detects unbroken hidden Z 2 x Z 2 symmetry, given 
as 



(On) 
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lim 

-3\-K 



exp 




(12) 



In fact, this was introduced in Ref. [SJ as a "parity cor- 
relation function" and shown to vanish in the Haldane 



phase but non- vanishing in a trivial phase. Here we dis- 
cuss Eq. (12) from a different viewpoint from that in 
Ref. M 

The non-local transformation (17 1 maps the non-local 
disorder parameter Eq. ( 12 ) to itself : 



(VOdV- 1 ) = lim /exp ji7r^e ,7rl: ^<'^Jf ) \ 

= ,4^( exp ( i7r § J '))' (13) 

where we used the fact that exp(±i7rJ 2 ) = exp(i7rJ z ) 
because J z only takes integer values. 

This correlation function can be interpreted as follows. 
The global 7r-rotation of spins (in the transformed basis) 
about z axis is a generator of the Z 2 x Z 2 symmetry 
Let us consider a localized operation, namely 7r-rotation 
about z axis only on the spins in the finite section bet- 
ween i and j. This "localized symmetry generator" is no 
longer a symmetry generator of the system. We apply this 
operation to the groundstate, and the overlap with the 
groundstate is measured. The limit \i — j\ —$■ 00 is taken 
afterwards. If the Z 2 x Z 2 symmetry is spontaneously 
broken, the application of the "localized symmetry ge- 
nerator" flips the order parameter in the finite section. 
Thus the overlap with the groundstate asymptotically va- 
nishes in the limit \i — j\ — > 00. Therefore, the disorder 



parameter (12) vanishes if the hidden Z 2 x Z 2 symme- 



try is spontaneously broken. This is quite analogous to 
the well-known disorder parameter in the quantum trans- 
verse Ising chain. 25 On the other hand, it does not vanish 
generically in a trivial phase where the hidden Z 2 x Z 2 
symmetry is unbroken. 



The discussion here implies that Eq. (12) acts as a 



disorder parameter for the hidden Z 2 x Z 2 symmetry, 
when the hidden Z 2 x Z 2 symmetry is well-defined. It 
would be the case even if the inversion (parity) symmetry 
is explicitly broken in the Hamiltonian, when the original 
argument in Ref. [Ml does not apply. (Although here we 
discussed the case of the tube, the same argument about 
the disorder operator applies to integer spin chains.) 



B. Edge states 

Possible quantum phases of the spin tube may be cha- 
racterized by the hidden Z 2 x Z 2 symmetry breaking (or 
non- breaking) . As in the case of single chain, sponta- 
neous breaking of the hidden Z 2 x Z 2 symmetry implies 
4-fold degeneracy of the groundstates, but only in the 
open boundary conditions. This implies the existence of 
the edge state (with half-integer spin, if the edge spin 
quantum number is well-defined.) 

It also implies that, we can investigate whether the hid- 
den Z 2 x Z 2 symmetry is spontaneously broken or not, 
by analyzing the existence of the edge states. If there are 
no edge states, the hidden Z 2 x Z 2 symmetry cannot be 



spontaneously broken. Existence of the edge state would 
suggest spontaneous breaking of the hidden Z2 x Z2 sym- 
metry. However, it should be noted that the edge states 
could appear by a different mechanism unrelated to the 
hidden Z2 x Z2 symmetry. 

The existence of the edge states can be analyzed easily 
in the strong-coupling limit (Jn <C J±). In the strong 
coupling limit, we can project to the groundstates of each 
rung, which changes according to the sequence Q , as the 
tube anisotropy parameter a is varied. 

Let us first discuss the S = 1 tube. In the isotro- 
pic regime (J11 <C J± and a ~ 1), each triangle tends 
to form singlets, and we thus expect a unique ground- 
state (with no boundary degeneracy) corresponding to 
the phase with unbroken Z 2 x Z 2 . As this phase has no 
spin at the boundary, it will be referred as the Sb = 
phase. On the other hand, in the anisotropic, "unfrustra- 
ted" regime (J|| -C J± and a < 0.5), the three spins of 
each triangle couple to form a J = 1 spin object. The re- 
sulting physics is essentially that of the spin-1 chain and 
we expect a ground-state degeneracy due to the boun- 
dary spins. In the language of the transformed Hamilto- 
nian this corresponds to the broken Z 2 x Z 2 phase with 
(O a ) ^ 0. We will call this phase the S b = 1/2 phase in 
relation with the spin 1/2 edge state of a spin-1 chain. 

We can also discuss the edge states in the weak cou- 
pling limit (J11 ^> J_l) with an heuristic argument. Taken 
separately, each chain of the tube is gapped, thus we can 
assume that a weak interchain coupling will not change 
qualitatively the physics in the bulk. On the contrary, 
solitary edge excitations of each of the three chains are 
expected to be very sensitive to any perturbation. As 
soon as a coupling Jj_ is introduced, they will be boun- 
ded, leaving a unique spin-1/2 degree of freedom at each 
edge. This again corresponds to the Sb = 1/2 phase, sug- 
gesting broken Z2 x Z2 . This picture is valid for basically 
any non-zero value of a. The only special point is a = 1 
where the translation symmetry in the transverse direc- 
tion leads to a bigger \ + \ degeneracy space as it happens 
with three spin-1/2 with identical AF couplings. 



Hidden Z2 x Z2 symmetry breaking in the weak 
coupling limit 



The edge state analysis implies that the hidden Z2 x Z2 
symmetry is spontaneously broken in the S = 1 tube, in 
the weak coupling limit J± <C J\\ for any value of a. 

However, there is a subtle issue in the weak coupling li- 
mit. It can be shown that the string order parameters ^ 
and (pi) exactly vanish at the decoupling point J± = 0. 
At this point, the groundstate is given by a simple pro- 
duct of the groundstates of each chain. The string order 



parameter p| can be decomposed as 
(O z ) = Jim ( S? a exp linJ^SfA S' kA 
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where () a is the expectation value with respect to the 
groundstate of chain a. While there are 3x3 = 9 terms, 
each one of them contains at least one factor of 
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2_^ S l,a 
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/= 



This is nothing but the disorder operator for the single 
chain, introduced in Ref. [2H and discussed in Sec. |IH A| 
It vanishes because the ground state of each chain in the 
Haldanc phase. As a consequence, the string order para- 
meter (18| for the tube also vanishes, apparently implying 
that the hidden Z 2 x Z 2 symmetry is unbroken. 



On the other hand, the disorder operator (12) also va- 



nishes because it is simply a product of the disorder ope- 
rators for the three chains. This rather suggests that the 
hidden Z2 x Z2 symmetry is spontaneously broken, in 
agreement with the edge state analysis. 

The resolution of this apparent contradiction is as fol- 
lows. The hidden Z 2 x Z 2 symmetry is indeed broken 
spontaneously even at the decoupling point J± = 0. Ho- 
wever, the string order parameters fl8| and ([9]), which are 
transformed to the ferromagnetic order in J z ' x by the 
non-local unitary transformation tffh, are not "good" or- 
der parameters to detect the symmetry breaking near the 
decoupling point. 

In order to detect the hidden Z2 x Z2 symmery breaking 
around the decoupling point, we introduce the following 
variation of the string order parameter : 

(O zzz ) = n lim (shSZzS^tt J ' S^S^Si- 

\k— 1\— >-oo \ 

The nonlocal transformation ftf\ maps this order para- 
meter to a local order parameter 



(y zzzy- 
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This does measure spontaneous breaking of the Z2 x Z2 
symmetry because S^S^S^ is odd under the global 
7r-rotation about x axis. 

At the decoupled point J± = 0, Eq. ( 16 1 reduces to the 



product of the standard string order parameter |6]) for 
the independent chains. Since Eq. ([6]) is non-vanishing in 
the groundstate of each chain, the "product" string order 
parameter ( 16 ) is also non- vanishing in the tube at the 



decoupling point. Therefore, the hidden Z2 x Z2 symme- 
try is indeed spontaneously broken even at the decoupled 
point Ji = 0, although the string order parameters $h 



and Q cannot detect the symmetry breaking. The edge 
spin Sb = 1/2 in the weak coupling limit, as well as in 
the strong coupling limit with anisotropy a < 0.5, can 
be understood as a consequence of the hidden Z 2 x Z 2 
symmetry breaking. Thus no phase transition is expected 
between these two regions, as both of them would belong 
to the hidden Z 2 x Z 2 symmetry broken phase. 

We note that the "product" string order parame- 
ter (16) is very similar to the O4 defined in Eq. (19) 
of Ref. 10. However, there is a crucial difference between 
the two-leg 5 = 1 ladder case studied in Ref. [10] and 
the three-leg ladder/tube case discussed in this paper. In 
the present case, the "product" string order parameter 
detects spontaneous breaking of the hidden Z 2 x Z 2 sym- 
metry, thanks to the relation eq. (17). However, in the 
two-leg ladder case, we find : 



VOaV 



Oa. 



(18) 



Thus 04 introduced in Ref. 10. does not detect hidden 
Z 2 x Z 2 symmetry breaking. 

In general, the product of string order parameters of 
each chain is an order parameter for the hidden Z 2 x 
Z 2 symmetry breaking in integer-spin ladder/tube with 
an odd number of legs, but not with an even number of 
legs. As a consequence, the hidden Z 2 x Z 2 symmetry 
(defined with respect to Eq. Q) is broken in the weak 
rung coupling limit of the odd-leg ladder/tube with an 
odd integer spin, but remains unbroken in the same limit 
if either the spin or the number of legs is even. 



D. Conjectured phase diagrams for S = 1 and 5 = 2 




1/2 1 3/2 2 5/2 



Figure 4: (color online) Numerical phase diagr am for the spin 
1 tube obtained with DMRG (see part VI B I. Both phases 
can be distinguished according to the string order parameters 
O z (see Eq.Q or |l6|) and part |VIBl[ ). The label S b for 
the different phases stands for the value of the spin of the 
boundary state for OBC. The transition between the Sb = 
and Sb — 1/2 regime is analyzed in more detail in this paper 
and turns out to be first order. Phase boundaries obtained 
from an effective model (see part |VI A[ ) are also shown with 
dashed red lines and agree quite well with numerical results. 



Let us now discuss the phase diagram. Here we propose 
the simplest phase diagram consistent with our analyses 
in the previous subsections. In Figure [4j we show the 
conjectured/numerical phase diagram for a S = 1 spin 
tube as a function of J11 and a. In one of the two phases, 
the hidden Z 2 x Z 2 symmetry is spontaneously broken 
and the edge state with Sf, — 1/2 appears. Although 
there are more edge state degeneracy on the special lines 
a = 1 and J± = 0, these lines are also a part of the 
broken hidden Z 2 x Z 2 symmetry phase. 

This phase diagram is also confirmed by numerical si- 
mulations ; the details will be given in Sec. |VIB| 
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Figure 5: (color online) Sketch of the phase diagram for the 
S — 2 spin tube. Again, the label Sb for the different phases 
stands for the value of the spin of the boundary state. The 
Sb — 1/2 phases (light red) must be separated from the others 
with real phase transitions (full lines). On the other hand, 
there is no reason to expect a real phase transition between 
the regions Sb — and Sb — 2 so we have plotted these 
boundaries with dash-dotted lines. 



We can also conjecture the phase diagram for the 5 = 2 
spin tube, as in figure [5] Based on the strong coupling 
analysis, we in principle expect three different phases : a 
singlet phase with no edge states, centered around a = 1, 
one phase with 1/2 boundary states, and two phases with 
Sb = 1 boundary spins. We can also discuss the weak 
rung coupling limit J± — >• in terms of edge states, by 
repeating the calculation of the decoupled triangle of the 
preceding section but by this time reasoning on the edge 
5 = 1 spins of the three 5 = 2 chains. We conclude, 
at the decoupling point Jj_ = 0, that there are no edge 
states for 0.5 < a < 2 and that there should be a single 
spin Sb = 1 at the boundaries for any other value of the 
anisotropy parameter. 

The edge state with 5& = 1/2 would correspond to a 
spontaneous symmetry breaking of the hidden Z 2 x Z 2 
symmetry, which clearly characterizes a distinct phase 
from those with Sb — 0,1. The phase with 5& = 1 also 
appears to be different from the Sb = 0, concerning the 
edge state. However, there is in principle no way to dis- 
tinguish the phases with 5f, = and 5^ = 1. This is 
suggested by the fact that hidden Z 2 x Z 2 symmetry is 



unbroken in the 5 = 2 Haldane "phase". 23 In fact, it 
was pointed out recently in Ref. [22] that the 5 = 2 Hal- 
dane "phase" with S^ = 1 is adiabatically connected to 
a trivial phase with Sb = 0. 

In our problem, the 5& = and Sj, = 1 "phases" are 
certainly adiabatically connected at the decoupled point 
J± = 0, where the system is just a collection of three 
5 = 2 chains with the Haldane gap. With an infinitesimal 
coupling J±, the gap should not close for any value of a. 
Therefore we expect that 5j> = and Sb — 1 "phases" 
actually belong to a single phase in which the hidden 
Z2 x Z2 symmetry is unbroken. 





Figure 6: (color online) The low energy modes deduced from 
the linear spin wave theory are connected to the three rota- 
tions of the initial triad around the axis x, y, z. When propa- 
gating to the tube, they underpin slow twists of the original 
structure which asymptotically cost no energy at long wave- 
lengths. 



IV. THE LARGE S LIMIT : NON LINEAR a 
MODEL 



A. Long wavelength description of the spin tube 

Having first examined the system from the strong cou- 
pling perspective, we now shift to the examination of the 
large-5 approaches whose greatest achievements culmi- 
nate with the non linear sigma model. The latter has 
shown to be particularly important in order to distin- 
guish the nature of the ground state, the low energy ex- 
citations and the possible critical points of an antifer- 
romagnelpSn23. It thus proves valuable to conduct such a 
study here to complete our previous analysis. The NLctM 
can be derived from the Heisenberg model when the spin 
S is large. In principle, it does not make any distinc- 
tion between integer and half-integer spins, as S is just 
one between multiple parameters allowed to flow conti- 
nuously to their renormalized values at long wavelength. 
However, the parity of the spin profoundly influences the 
value of the Berry phase, a purely quantum quantity ori- 
ginating from non zero overlaps between coherent states 
and entering the NLctM action. As shown by Haldane 1 , 
the value of the Berry phase eventually governs the pro- 
perties of the system in the infrared limit. 

Although the spin-wave expansion cannot make good 
quantitative predictions on a magnetic model in one di- 
mension, it is useful to carry this analysis in order to 
identify the low-energy, long-wavelength degrees of free- 
dom in the spin tube. These degrees of freedom will help 
later on to construct a well-defined order parameter for 
the NLtrM. A standard spin wave analysis shows that 
there are three low energy modes depicted in Figure ([6]) 
and their canonical conjugate. With the expression of 
the low-energy modes e? , e< and e< , and their respec- 
tive conjugate modes e¥, ef and e?, we can reexpress the 
slowly varying spin degrees of freedom. The spin opera- 
tors can be rewritten in a compact form by introducing 

an infinitesimal SO (3) rotation matrix TZi — exp I irriiJ ) 
with J standing for the generator of 50(3), and a vector 



Lj if we identify 
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The original spin operators read very simply : 
S i>a = (-lynma + (Li - (U • n Q )n Q ). 



(19) 



where the vectors n a are given in |2]). The number of 
degrees of freedom on the left and on the right hand side 
of this equation is the same. Thus, this operation can be 
regarded as a simple change of variables, but if one is 
interested in a long wavelength limit, the new variables 
TZi and Lj are the most adapted to describe the physics 
of the anisotropic spin tube. 



B. Derivation of the NLctM 

After having derived the low energy modes of the 
theory, we can now focus on the construction of the 
NLsM. The first point is to define a suitable local or- 
der parameter. In the case of a coplanar configuration, 
the role of the order parameter will be played by the 
50(3) rotation matrix. The underlying idea behind the 
construction of the NLsM is then to consider the system 
in its symmetry-broken phase and to take into account 
quantum fluctuations around the direction of the order 
parameter. It is actually not necessary that the system 
admits a broken-symmetry phase (it has none in Id), it 
just needs to be at least ordered locally, which is the case 
in the weak rung coupling limit J11 ^> J[_. The conti- 
nuum limit is then reached within a Hamiltonian for- 
malism or a path integral (lagrangian) approach 29 . One 
of the greatest possibilities given by the NLsM is that 
it can be investigated with renormalization group tech- 
niques for which one can determine the nature of the 
possible fixed points governing the physics at long wa- 
velengths. The quantum NLsM for triangular geometries 



and its RG anal ysis h ave already been extensively studied 
by Azaria et a/j 18 * 30 l We shall refer later to their work 
for the characterization of the spectrum. For now, our 
starting point will be different. We will build the NLsM 
in the lagrangian formalism, following the construction 
of Dombre and ReacP-^. This path integral approach is 
particularly illustrative regarding the construction of the 
Berry phases. These complex terms will play a major role 
in our analysis of the quantum spin tube. 

In the long wavelength limit, the fine scale of the lat- 
tice becomes irrelevant and the lattice spacing A can be 
taken to zero. To obtain a well-defined continuum limit, 
we need to choose a local order-parameter field which has 
a smooth spatial variation on the scale of A. To do so, we 
make use of the following ansatz^ : 



side, 



S' BP = iSj2(- 1 ) i "lftina} 



(24) 



z,a 



takes a particular significance when one allows for the 
possibility of singularities in the action. These singulari- 
ties are naturally present in the system because we start 
from a lattice description with discrete variables and not 
fields. However, we wish not to take them into account 
now and we will let apart this term for the moment. We 
will reconsider it when evaluating the role of the topolo- 
gical defects in the theory. 



Si,„(t) = S 



Ki(t) ((-l) 4 n Q + ALi(t)) 



^(*)((-l) 



AL*(*)) 



(20) 



where the fields H and L depend on the time t and the 
lattice coordinate i. This is identical to ( 19 ) at first order 



The second member of Eq. ((22 1) is nothing but the 
Hamiltonian (fTl) where the quantum spin operators have 
been replaced by the ansatz (21 ). The Hamiltonian part, 
at second order in A, reads : 



Sh = I dxdT ( -y Striae! - 6 AS" TL(x, t) ■ L(x, t) 



in A. Here however, to make a coherent calculation, we 
will need to keep the development up to second order in 
A : 



+ S Tr[P(K{x, T)~ l d x Tl{x, t)) 2 



with 



S,, a « ^. i ((-l) t n a + A(L l -(n Q .L l )n a ) 
// T 2 9 

+ {-lyx 2 ' 1 



(21) 



Str. 



iad 



o Z^ S *.a 



1-a 



(Si.i 



(n .Li) n Q - (n .Li)L, 



In particular, the square of the magnetization of the 
whole triad is given by : 



eM -* 2 



a — 1 



9A 2 S 2 (TLi 



X 2 S 2 



a- 1 



1 — a *■ — ' 



where T a p = 5 a/3 - lJ2a n aan a p, and (TL) a = 
S/3 TupLp. We set again Jy = J± = 1 for simplicity and 
concentrate on the effects of the anisotropy parameter a. 
The action we wish to estimate iipS : 



•Snlo-m = Sbp + Sh 

= i 



X;«[S 4>0 ]- JdrH[S iia ] 



(22) 



Since the action is quadratic in the field L(x,t), we 
can integrate the field out and finally express the action 
solely in terms of the SO (3) matrix field lZ(x,t) : 

Snl*m = S fdxdr {Tx[P{Tl{x,T)- 1 d x '1l{x,T)f] 

+ ll:[Q(^(x,r)- :L a r ^(a:,r)) a ]) + iS' BP . (25) 

Here P and Q are diagonal matrices whose expression is 
better given by the spin-stiffness and susceptibility ten- 
sors^ : 

Xa/3 = -Tr (Qt a tp), p a = -Tr{Pt a tp), (26) 

With Xaf3 = XaSa/3, p a = PaSa/3 and I 

S(l + 2a 2 ] 2 ^ _ 95a 

' A(2 + 8a(4 + a))' 



t denoting imaginary time. The first part of the action is 
the Berry phase term. It measures the total area covered 
by each of the spins Sj, a on the sphere of radius S. Up 
to second order, the Berry phase term reads : 



Xi = 
X3 = 

Pl = : AS' j 1 



Aa(l + 4a(3 + a + 4a 2 )) 
Sfl_ 8 

A \a -1 + 4a(2 + a) 
1 
2a^ 



X2 



P2 = 3XS, P3 = \S[2-^ )(27) 



1 
2a 5 



i.a 



(-l)^[^n a ] + i3S I drdx (TU) ■ V, 

(23) 
where V a = -\e a ^{1Z- 1 d T 1l)p 1 and e a ^ is the totally 
antisymmetric tensor. The first member of the right-hand 



Here t a are the generators of the SO (3) group. Finally, 
we introduce the fields 7Z(x,T)~ 1 d IJ ,lZ(x,T) — oj^t a . The 
action reada^ : 



5* 



S 



NLo-M 



drdx 



( a B 

[X a f3U W 



• p a pUJ x 0JZ 



(28) 
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C. Bare analysis of the NLuM. 



(|25| and (|28j), 
rould be able 



The two formulations of the action, eqs. 
are valid for all a > 0.5. In particular, we si 
to recover the isotropic limit a = 1, and the unfrustrated 
cases a — 0.5 and a — » oo. For a = 1, one notes that 
Pi = pa 7^ P2 and Xi = X3 7^ X2- In the language of the 
50(3) matrices, this translates into an additional 5*0(2) 
global right symmetry of the action : 1Z — > TWr. This 
symmetry is reminiscent of the discrete C% v symmetry 
of the triangle for a = 1. Since the configurations of 
fields classically minimizing the action also possesses a 
50(2) symmetry, this model is referred to as the 50(3) x 
50(2)/50(2) NLo-]VPnSI. When a = 0.5, one finds Pl = 
Pi-, Xi = Xi an d X3 — P3 = and recovers the description 
of the collinear antiferromagnet in terms of a 0(3)/0(2) 
NLcrM. 

There is another, third representation of the action 



(25)- (28 1 that nicely illustrates the effect of the aniso- 



tropy parameter a. Remembering that a 50(3) matrix 
is nothing more but a set of three orthonormal vectors : 
(e&)a = TZ a b, we can use the fact that e 2 = e,3 x ei to 
rewrite the bare action in terms of two orthonormal unit 
vectors : 



5, 



5nlctM 
5i 

coupling 



5i + 53 + 5 COU pij llg + S BP , (29) 

— / dxdr (c a {d T e a ) 2 + ^{d x e a ) 2 j , 



2fl. 



da;dr(e 3 • d T ei) z 



where the new constants can be easily expressed as a 
function of the spin-stiffness and susceptibility tensors. 
The behavior of the different couplings as a function of 
a are easily obtained from the expression of the spin- 
stiffness and susceptibility tensors. An important point 
is that for a — > 0.5, g\ — > 00 and k —¥ i.e. when a = 0.5 
the ei field becomes a spurious degree of freedom with 
null stiffness ! This is consistent with the collinear picture 
in the range < a < 0.5. Such a model is in fact well 
described by the fluctuations of a single unit vector 63. 



D. Renormalization group analysis of the NLcrM. 

The above bare analysis of the NLcrM is insufficient 
to describe properly the behavior of the spin tube in the 
infrared limit. As it is well known from the study of the 
quantum spin chain, quantum fluctuations always renor- 
malize the parameters entering a NLcrM action and even- 
tually drive the system into a quantum disordered state 
in lcp2l. Thus, in order to understand the properties of the 
system at long wavelength, one must perform a renorma- 
lization group analysis to determine how do the different 



coupling constants (27) renormalize. We aregoing to use 
the results for the one-loop RG equation d 18 * 30 *, starting 
with the set of bare couplings (27). To make the dis- 



troduce the two anisotropy parameters a 2 = 1 — Pi I ' Pi 
and 03 = 1 — pa/pi and the coupling g = 2/ pi, the latter 
playing the role of an effective coupling constant. The set 
of couplings 7 = {01,02,03,02, 0:3,(7} obeys the general 
RG equations : 



In 



= ~/3({7}) 



(30) 



We have integrated numerically these equations for dif- 
ferent values of a ranging from a = 0.55 to a = 0.95 
(similar behaviors were also observed for values above 
a = 1). The numerical integration of the RG equations 
yields the unambiguous result that the symmetry is dyna- 
mically enlarged in the infrared limit, similarly to the hi- 
gher dimensional cases. For any value of a > 0.5, the spin 
wave velocities renormalize to the same value c* — c^ — 
c 3 while the two anisotropy parameters a-i and 03 fall to 
zero. The symmetry of the model in the long- wavelength 
limit is therefore 50(3) x 50(3)/50(3) sa 50(4)/50(3). 
The coupling g diverges, as one could have expected in 
one dimension. Hence, there seems to be no qualitative 
differences between the cases a = 1 and a > 0.5, a 7^ 1 
at the one loop level, suggesting that a deviation from 
the point a = 1 is an irrelevant perturbation. However, 
we have not taken into account so far the role of the 
Berry phase term, which as we are going to see plays an 
important role when a^l. 



V. BERRY PHASES AND INSTANTONS 

A. Instantons in SO (3) NLcrM 

The continuous part of the sigma model does not make 
any distinction between the integer and the half-integer 
quantum spin tube. In fact, the preceding RG equations 
suggest that the model is gapped in both cases, and ad- 
mits a unique ground state for any a. Nonetheless, the 
DMRG data show unambiguously a dimerization of the 
ground state of the spin tube for 5 = \ at a = l | 14 | 33 l 
Analogously, the Majumdar-Gosh model, whose NLcrM 
also has t he 50(3) x 50(2)/50(2) symmetry, is dime- 
rizecpSEll. In a single spin chain, the difference between 
integer and half-integer spins can be explained in the 
NLcrM by the presence of a topological term in the ac- 
tiorffl Here, such a term is absent because of the trivia- 
lity ofthe second homotopy group of the 50(3) mani- 
lolrpEl: 

tt 2 (5O(3))=0 

However, only continuous space-time configurations of 
the field lZ(x, r) have been considered up to now. In 
fact, there also exist configurations containing vortices 
with singular cores. These defects originate from the non 
trivial first homotopy group of 50(3) : 



tinction between isotropic and anisotropic cases, we in- 



7Ti(SO(3)) 



(31) 
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For classical antiferromagnets on the triangular lattice, 
these vortices are argued to be the driving force of a phase 
transition^. 

In quantum systems, topological defects radically 
affect the behavior of the disordered phases of the 
0(3)/0(2) NLerM in 2cP, leading the system to dimeri- 
zation, and of the 50(3) x 50(2)/50(2) NLcrM in lcP. 
The specificity of our system is that the 50(2) symme- 
try is, at least at the bare level, no longer present when 
a^l. Thus, we would like to investigate the conjugate 
action of the topological defects, also known as instan- 
tongS2l, and of the anisotropy in the spin tube. For integer 
5, we will see that the presence of the topological defects 
gives rise to the emergence of 25 peculiar values of a, 
that we could associate with the critical points determi- 
ned from the strong coupling approach. 

We would like to review first the nature of the ins- 
tantons in our system. Instantons are topological defects 
associated to the symmetry group of the order parameter. 
The 50(3) group manifold is isomorphic to a ball of ra- 
dius ir in three dimensions whose diametrically opposite 
points on the surface are identified. One can associate to 
a rotation around an axis n by an angle 9, the vector #n 
with 9 e [— 7r,7r]. The redundancy between two opposite 
points on the surface of the sphere stems from the iden- 
tification between a rotation about an axis n of angle it 
and the rotation about the same axis of angle — -k. It is 
then clear that the 50(3) manifold is non simply connec- 
ted, with the ensemble of closed path in 50(3) dividing 
into two classes : one containing the loops shrinkable to a 
point, the other ones containing strings joining two oppo- 
site points of the ball. This is equivalent to say that the 
first homotopy group of the 50(3) manifold is given by 
(31 1. Considering the evolution of a matrix IZ(x) through 
space, an element of the non trivial class is : 



cos9(x) sin9(x) 

H(t) = [ -sin 0(a) cos0(a:) 

1 



(32) 



with 9(x = 0) = and 9(x = L) = 2tt, where L is the size 
of the system. Conversely, the trivial class will consist of 
matrices which stay close to the identity matrix at all 
positions. 

Turning back to our 1 + 1 dimensional problem, sup- 
pose now we start from a configuration in the trivial sec- 
tor with all ei vectors pointing up and all e2 vectors 
pointing right, where again lZ a p = (e a )^ (see Fig [7]). If 
nothing "sudden" happens, i.e. if the time evolution pro- 
cess is sufficiently smooth, the chain should visit other 
configurations but stay in the trivial topological class. 
However, it is also possible that some non trivial configu- 
rations arise during time evolution that will connect the 
two classes of path. These are the instantons. A pair of 
instanton (going from the trivial to the non-trivial class) 
and anti-instanton (i.e the opposite) is represented on 
Fig. [7] In the continuum, an instanton will appear as a 
singularity. It is clear that such an event is unlikely to 
happen if the tube is ordered. However, since the model 
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Figure 7: (color online) Evolution of the spin tube from the 
two different topological sectors of 50(3). Each triad is repre- 
sented by two orthonormal vectors ei(x, t) and e2(x,t) that 
one can connect with the classical configuration of spins on 
each triangle (Si,i, Si,2, S;,3/^. In the continuum, ei and e2 
also stand for the first two vector columns of the rotation ma- 
trix lZ(x,t). Starting from a Neel configuration, the system 
tunnels to a non trivial configuration via an instanton event 
(+). The system returns to the trivial configuration via an 
anti-instanton (— ). 



is disordered at long wavelengths, these events will even- 
tually proliferate. Now, it may be that the proliferation of 
instantons is constrained by the Berry phase term. Here, 
we would like to calculate : 



S'rp = iS\^ u)[Hia. a ](- 



7 BP 






■iy 



(33) 



which is the discrete part of the total Berry phase ( 23 ) 
that we let apart. For this purpose, we follow Dombre 
and Read again and consider a first time path 1Z(t) sa- 
tisfying the closed boundary conditions and a second one 
W(t) = K(t)+5K(t) infinitesimally close to K{t). The 
difference of Berry phases between the two paths can be 
easily evaluated to be : 

SSgp = iS /dT^((5ftn a ) • {d T TZn a x TZn a ) 

J a 

= -iS f drv^n-^n)^ lj2n a _ 



P' 



B. Isotropic case, a — 1 

In the isotropic case, we have the important result 
that : 



SS' BP = 



(34) 



and any smooth change in the history of TZ(t) will not 
change the value of the Berry phase of the triad. Thus, 
this quantity can be used to index the two classes of 
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7Ti(50(3)), exactly like the hedgehog number classifies 
the configurations of the spins in two dimension : ! 1 * 29 !. Be- 
cause the quantity is a topological invariant, we just need 
to calculate it for one path representing each class. For 
the trivial class, we can take the identity matrix so that 
J2 a w [7^(£) n a] = [47r5]. For the non-trivial class, we can 
consider the rotation of the triad around an arbitrary 
axis. In this case, the Berry phase of the triad will be 
2irS [AirS]. So, the alternating sum (33 1 reads : 



S' t 



BP 



iSy^^2wqi(- 



1) 



(35) 



a A 



c^ 



Ci 



C^ 



I x l 



d 



2nS -2nS 2nS -2nS 2nS 



where qi = 0, 1 depending on which class the matrix IZi 
belongs to. Consequently, the total Berry phase will be 
or 2ttS depending on the number of non-trivial paths. If 
S is an integer, the Berry phase has no effect. But if S is 
a half-integer, we see that there are two different values 
for Sbp, defining two different vacua. 
To see the influence of the instantons on the system, we 
remember the arguments of Rao and SerPB. An instan- 
ton is a discontinuity in the Berry phase of two neighbo- 
ring triads. Because of closed boundary conditions in the 
partition function, an instanton necessarily comes with 
an anti-instanton. As we saw, the creation of such a pair 
links the two vacua labelled by Sbp = and Sbp = 2nS. 
The instantons are situated on the links of the lattice (as 
they live on the plaquettes of the lattice in the (2+l)d 
case). A pair of instanton, anti-instanton defines a string 
of a given size. If this size is even, the Berry phase of the 
string is ; if it is odd, it is 2ttS (Fig[8]). It is then easy to 
see that if S is half-integer, there will be destructive inter- 
ferences between paths with strings of different sizes. In 
particular, if we shift an instanton by one lattice site, we 
expect the dynamical contribution from the Hamiltonian 
to not change, but the Berry phase to change by 2nS. 
For instance, the two paths of Figure [8] will contribute in 
the partition function : 



(1 



- e 2i7rS )e SN 



(36) 



For half-integer spins, different instantons-anti- 
instantons contributions are compensated by destructive 
interferences. Thus, the two topological sectors q = 0, 1 
are non connected and we are left with two degenerate 
ground states labeled by the two elements of n 1 (SO(3)). 
As shown by Read and Sachdev in a large N analysis of 
the (2+l)d Heisenberg modeP^, this kind of destructive 
interferences between instantons leads to dimerization 
in disordered phases. This seems to be the case here : 
the spin 1/2 isotropic model is known to be dimerized 
by DMRG. 

On the other hand, integer spins allow instanton events 
to proliferate as all events come with the same phase. 
This makes that the two vacua are well connected. This 
"tunneling" between vacua lifts the degeneracy and the 
ground state is therefore unique. 



Figure 8: (color online) Two space time configurations diffe- 
ring by the shift of one instanton. The two crosses represent 
the positions of the singularities. A straight line represents a 
trivial path in SO (3) while a loop is a non-trivial path. The 
Berry phase associated with each loop is ±2^5* depending on 
the sublattice. 



C. Anisotropic case, a/1 

For a/ 1, the difference in Berry phase between two 
matrices 1Z and 7Z+57Z belonging to the same topological 
class is : 



%p> 



6S' BP = -iS /"drV^ft-MftW^iv 

J a 

J dr Se 3 ■ (d T e 3 x e 3 ). (37) 



,51-q 
'2 a 



In the anisotropic case, the Berry phase can no longer 
be used to classify topological classes. For example, the 
path contribution to the partition function from the two 
configurations drawn on figure IS] is now : 



Z 



I e 2i7rSL ( e »f ^sr Jo dx I dT S x e 3 -(d r e 3 xe 3 ) /gg\ 

fg +1 dx f dr d x e 3 -(d T e a xe 3 )\ Snwm t 



e i2nS e i§^ >'■ < 



Here, we have been careful to write the total Berry 
phase as a continuous integral over the string. By doing 
so, we made the approximation that the order parameter 
1Z(x, t) is sufficiently smooth so that the derivatives d x e% 
are well-defined. This development is valid if we stay in 
a given topological sector of SO(3). For a configuration 
with many instantons, we should separate the contribu- 
tions from the different topological sectors and write it 
as a sum of integrals : 



s' BP = iY,^s qi (-iy 



(39) 



c i p - 1 

2 a ^ 



Xi + l 



dx / dr d x e 3 ■ (d T e 3 x e 3 ), 
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where {x±, . . . , Xi} denotes the position of the P instan- 
tons. Note that without instantons (i.e if 1Z(x, t) is a 
smooth field everywhere), we can regroup all the integrals 
into a single one and this term identically vanishes : 

/+oo p 

dx / dr d x e 3 ■ (d T e 3 x e 3 ) (40) 

-OO J 

dxdr [9 T (ei • d x e 2 ) - d x (e x ■ d T e 2 )} = 0, 



( 39 1 into a single one : 



given the periodic boundary conditions we imposed. 

Let us finally recover some well-known result in the ex- 
treme limits a — 0.5 and a — > oo. In this case, the sym- 
metry of the order parameter reduces to 0(3)/0(2) = S 2 . 
There are no instantons in this case since tti(S 2 ) = 1. It 
is then straightforward to show that ( 39 ) reduces to : 



Qtotf 

>~>BP 



s 

= % 2 



/+oo p 

dx I dr d x e 3 -(d T e 3 xe 3 ) = 2m Sn n € Z 
-oo J 



We recall that having a non-trivial skyrmion number for 
a smooth space-time configuration of the vector field e3 
requires discontinuities on the field ei. However, as this 
last field gets zero stiffness for a — 0.5 and decouples we 
recover for the field e3 the form of the NLcM with the 
correct topological term of a single chain of spins S as we 
should. 



1. Half-integer spins 

Reiterating the argument that led us to the twofold de- 
generacy of the ground state for a = 1 and half-integer 
spins, we find with (38) that the different instantons- 
anti-instantons contributions do not cancel out anymore. 



The tunneling process between the two topological sec- 
tors is present and the topological degeneracy is lifted. 
Consequently, we can make use of an important result for 
spin chains, the Lieb-Schultz-Mattis theorem*™ sugges- 
ting that the system is in a gapless phase. This theorem 
states that spin Hamiltonians with local interactions and 
an half-integer spin per unit cell like (fTl) either support 
gapless excitations or have a ground-state degeneracy. 
Ruling out the possibility of a degeneracy here tends to 
the scenario of a critical behavior. This is indeed what 
appears in the the study of T. Sakai et al where, for 
S = 1/2, the DMRG data point at a preservation of the 
spin gap only in a narrow range around a = lQ3l. 



2. Integer spins 

An interesting application for integer spins is a pos- 
sible extension of Haldane's conjecture to the quantum 
spin tube. Reconsidering again the Berry phase term, we 
examine the possibility of rewriting the sum of integrals 



p-i 



2 a 



8=0 

Q 



Xi+l 



i- ^ / dx dr d x e 3 ■ (d T e 3 x e 3 ) 
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/+oo p 

dx / dr d x e 3 ■ (d T e 3 x e 3 ). (41) 
-oo J 



hand side of ( 



Because of (40), we saw that the integral on the right 



1 ) must vanish for any smooth configu- 



ration of the field lZ(x,t). However, it is possible that 
the field e 3 (x,t) is smooth but that 1Z(x,t) is not (see 
for instance fig [7] : the vectors ei and e 2 change sharply 
of direction where the instantons take place but e 3 re- 
mains constant). In this case, writing the Berry phase as 
a single integral is allowed and this integral will be dif- 
ferent from zero. However, we must emphasize that the 
identification (41 ) is not totally complete, since we elude 
all the instantons events where e 3 is discontinuous. Ho- 
wever, in the region a ~ 0.5 we saw that in the bare 
action the stiffness of the e! is very small. One can then 
suppose that in this limit the low energy configurations 
with non-trivial topological index are those where the e 3 
is smoothly varying and the necessary discontinuities are 
in the e! field configuration. Keeping this picture even 
for larger values of a, from (41) one then recognizes a 



topological term for the unit vector e 3 . But this time, it 
is multiplied by a factor (a — l)/a. For the 0(3)/0(2) 
NLerM, it is known that such a term would lead to a si- 
gnificant change in the spectrum of the model if it is an 
half-integer, in which case the NLcrM is gapless. Here, we 
find 25 particular values of a for which this happen : 



S 



s-( P - 



S' BP = in{2p + 1) (42) 



-S <p+ - < S, peZ, 

the last inequalities coming from the condition a > 0.5. 
For a — a p , the Berry phase reduces again to an odd 
multiple of tt. Finally, the full anisotropic sigma model 
at these points read : 

Snl.m = [dxdr (^r- (c a {d T e a ) 2 + ^(d x e a )A (43) 



+ n(e 3 ■ d T ei) 2 + i 



2p+l 



d x e 3 ■ (d T e 3 x e 3 ) 



Would the RG analysis of the preceding section have 
predicted a decoupling of the field ei, we would have 
concluded that the last equation represents 2S critical 
field theories, each corresponding to an SU (2)i Wess Zu- 
mino Novikov Witten (WZNW) model, as it happens for 
a dimerized spin S chain 39 . However, the RG results sug- 
gest the opposite scenario : the coupling between the ei 
and the e 3 is a relevant perturbation and any non-zero 
skyrmion configuration must come with a fugacity, cor- 
responding to the energy cost of a discontinuity of the 
ei field configuration. This would exclude the scenario 
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for a SU(2)i WZNW criticality. Although the one-loop 
RG results are well suited to study the vicinity of the 
point a = 1, one can question their validity in the quasi- 
collinear regime a ~ 0.5. The nature of the transition 
points J — > J + 1 (J is the total spin per triangle, see 
equation Q is thus unclear for the case J big (i.e. close 
to a — 0.5) and we are going now to study the transition 
J = to J = 1 in the case of an S = 1 tube. 



VI. 



THE S = 1 CASE 



In this section, we will focus on the spin-1 case for 
which two quantum phase transitions are expected (close 
to a = 1/2 and a = 2 respectively when J\\/J± ^C 1). 



A. Effective model for the S — 1 tube 

In the small Jy limit, we can apply simple perturbation 
theory in order to obtain an effective model which should 
be valid close enough to a critical point. As recalled in 
Sec. |TlJ a single triangle exhibits at low-energy a level 
crossing between one singlet and one triplet states both 
at a = 1/2 and a = 2. If we restrict ourselves to the 
neighboring of one level crossing, then we can build an 
effective model by keeping only these low-lying degrees of 
freedom. Because the Hilbert space of one singlet plus one 
triplet is equivalent to two spin 1/2, we prefer to describe 
the effective model in terms of effective spin 1/2 variables 
so that the effective model of the spin-tube becomes a 
spin-1/2 2-leg ladder hamiltonian. 

By performing first order (in J» ) degenerate perturba- 
tion, we end up with a SU(2) spin-1/2 ladder that only 
contains 2-spin exchange interactions of the form : 

^eff = Z_^, J-iSiA ' ^.2 + J\\(Si t i ■ <Si+i,i + Si+i,2 ■ Si+1.2) 

i 

+ Jd'$2(Si,i'Si + i,2 + Si + i,2'Si,i) ( 44 ) 

i 

and the effective exchange are as follows : 

a -1/2: J ± =2a-1, J {] = (~ + ^)J h 



Jd = (y " ^|| ~ 0.82J|, 



a ~ 2 : J ± = 2 - a, J11 



Jd = 77 -Til ~ 0.56 J,, 
y 
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J\\ . 



(45) 



(46) 



This mapping allows for a straightforward explanation of 
the occurence of quantum phase transitions. Varying a is 
equivalent to changing the effective rung exchange from 
strongly positive to strongly negative, which means that 
the spin-1/2 ladder is in a rung-singlet phase on one side 
and in a Haldane phase on the other side. Because this 
spin-1/2 model is simpler and has already been studied 



intensively, we can use some results from the literature 
to clarify the nature of the phase transition. 

From the bosonisation point of view, which is valid 
when Jy is the dominant energy scale, Nersesyan and Ts- 
velik have argued that there should be a transition when 
J± = 2Jd with the possibility of deconfined spinons-^! A 
more refined analysis by Starykh and Balents has shown 
that marginal interactions modify these conclusions so 
that the transition between rung singlet and Haldane 
phase becomes either first order or has an intermediate 
columnar dimer phased These estimates J± = 2Jd for 
the quantum phase transition are plotted on the phase 
diagram in Fig. [4] 

From the numerical point of view, early DMRG simu- 
lations^ were interpreted in favor of a second order (res- 
pectively first order) phase transition for Jd/J\\ smaller 
(respectively larger) than 0.287. The absence of an inter- 
mediate dimerized phase was confirmed by more recent 
numerical work^ 2 although these studies do not agree 
on the order of the transition : either it is always first 
ordeiH'or it could be continuous for small Jd/ J\\ ^ 0.2=2! 
Given that our effective models have a relatively large 
ratio Jd/J\\ (respectively close to 0.42 and 0.38 for both 
critical cases a ~ 1/2, 2), all numerical studies agree that 
the phase transitions are first order. 



B. DMRG results for the S = 1 tube 

In order to have an unbiased answer, we have decided 
to perform numerical simulations of the 5 = 1 tube with 
the powerful DMRG algorithm^ for several values of Jm. 
Simulations are done mostly with open boundary condi- 
tions (OBC) with system sizes up to 3 x 64, but also with 
periodic boundary conditions (PBC) on some cases. Ty- 
pically we keep up to 1600 states, which is sufficient to 
have a discarded weight smaller than 10~ 8 . 



1. String order parameter 

In order to draw a numerical phase diagram and to 
compare it with the conjectured one (see Fig.Rl), we have 
computed the z-component of the string order parameter 
(see Eq. (l9|) for several values of Jy and a. In order to 
extract the bulk value and avoid finite-size effects due to 
the edges, we have taken the following definition in our 
simulations : 



(O z 



Jl/i ex P 




3L/4+1 



(47) 



In Fig. [9j we plot this quantity as a function of the 
frustration a for a small J\\/J± =0.1. We conclude that 
the string order is finite for a < 0.57 and a > 1.92, and 
it vanishes elsewhere, i.e. our model does exhibit quan- 
tum phase transitions. Therefore, as a function of a, we 
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Figure 9: (color online) String order parameter along z of 
the spin-1 tube as a function of the anisotropy parameter for 
J\\ — O.lJx and various lengths L. The extrapolated spin gap 



(see Fig. 13 below) is also shown, as well as the non-local 
order parameter of the spin-1 Haldane chain. 
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Figure 10: (color online) String order parameter along z of 
the spin-1 tube as a function of the anisotropy parameter for 
Jii = J± and various lengths L. A linear extrapolation with 
the largest sizes is also plotted and is finite for all a. 



find successively a topological phase (with Sb = 1/2 in 
the presence of OBC), a non-topological one, and again 
a topological phase (with Sb = 1/2 in the presence of 
OBC). This is the behaviour expected from the pertur- 
bation and the mapping to an effective spin-1 or spin-0 
chain. Indeed, for small or large a, we can derive an ef- 
fective spin-1 Haldane model for which the string order 
parameter is knowrP^to be ~ 0.374, which is close to our 
value in both limits. 

In order to complete our phase diagram in Fig. |4j we 
also compute the string order parameter for larger J\\/J± 
where perturbation is no more valid. Data are shown in 
Fig. [lO] and have a quite different behaviour : now, the 
string order parameter is finite for all a, i.e. we have no 
phase transitions along this line. 



0.06 



0.04 



O z 



0.02 




Figure 11: (color online) String order parameter along z of 
the spin-1 tube as a function of J\\/J± for a fixed anisotropy 
parameter a — 0.75 and various sizes. 



By computing (O z ) for various J» and a, we estimate 
that the tip of the J = lobe occurs for J\\/J± = 0.67 
and a = 1. 

Now, we present data for a vertical cut in the phase 
diagram of Fig. [4lby fixing a = 0.75. By varying J\\/J±, 
the string order plotted in Fig. fill vanishes for J\\/J± < 
0.34 and is finite beyond. In Fig. [TTJ the string order 
parameter (O z ) remains finite up to J\\/ J±_ ~ 3. However, 

" (O z 



IIIC 



vanishes in the 



as we have discussed in Sec. 
weak coupling limit. 

In order to ascertain that there is no other phase 
transition when going to the decoupled chain limit, we 
have computed the product string order parameter of 
Eq. ( 16 ) as well as the usual one for the non-frustrated 



case a = 0. Data are shown in Fig. [12] for both order 
parameters. When the chains are almost decoupled, the 
product string order parameter is close to the product of 
the standard string order parameter for the independent 
chains. On the opposite side, when the dominant cou- 
pling is J_l, the pe rturba tive argument that we have gi- 
ven above (see Sec. VI A | indicates that the 3-leg ladder 



behaves effectively as a spin-1 chain for which the string 
order parameter (O z ) reduces to the usual one; this is 
indeed what is found numerically on large system size. 



2. Spectral gap 

We have also calculated the excitation gap by DMRG. 
In order to estimate the bulk gap, the gap is extrac- 
ted differently from the finite-size spectrum depending 
on the boundary conditions. For open boundary condi- 
tions (OBC), in a Haldane- like phase, the real gap should 
be calculated between the S = 1 sector and the S = 2 
sector since the sectors 5 = and 5 = 1 are already 
degenerate because of the edge states. On the contrary, 
for a singlet-like state, the gap is defined between the 
5 = sector and the 5=1 sector. For periodic boun- 
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Figure 12: (color online) Various string order parameters of 
the spin-1 ladder (a = 0) as a function of J±/J\\ for va- 
rious sizes. Usual (respectively product) string order along 
z is shown with open (respectively filled) symbols. The da- 
shed line indicates the known value for the spin-1 Haldane 
chairP. 



dary conditions, there are no edge states and thus, the 
gap is uniquely defined to be between S = and 5=1. 



The evolution of the gap is presented on Figure 13 (left 
for OBC and right for PBC). For OBC on largest system 
sizes and fixed J\\/J± = 0.1, the DMRG indicates that 
the gap between 5 = and S = 1 is almost zero for 
a < 0.57 or a > 1.92 but is finite in between (data not 
shown). The gap between 5=1 and 5 = 2 is plotted 
~ for Jii 



on Fig 
rence in 
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O.lJi and exhibits a striking diffc- 
tnree regions. For a < 0.57 or a > 1.92, the 
gap is finite and almost constant with a. In contrast, in 
the intermediate a region, the gap increases almost li- 
nearly away from these critical points. Extrapolation of 
the data for large system sizes seems to go in favor of a 
finite gap everywhere (see Fig. 13), but note that the ex- 
trapolated gap at the critical points is extremely small. 
Both critical points are identical to the values we had 
found with the string order parameter. We observe that 
the gap is roughly constant in the Sb = 1/2 phase, except 
in the vicinity of the transition point. Note that this is 
in accordance with the qualitative picture of the strong 
coupling limit. For small or very large a, the effective 
model is a spin-1 chain with an effective spin exchange of 
order J» (but independent of a) ; therefore, the spin gap 
essentially depends on the value of Jy but is independent 
of the anisotropy parameter a. 

In order to make connection with the perturbation 
theory, we have also performed simulations of the effec- 



tive spin- 1/2 ladder (see Eq. (44)) corresponding to our 



parameters choice. Due to the Hilbert space reduction, 
we are able to simulate larger clusters. As can be seen on 
Fig. 14 we obtain a very good agreement between both 



sets of data since we are indeed considering a small J\\ 
case where perturbation is expected to be accurate. We 
have plotted an infinite-size extrapolation by using the 




Figure 13: (color online) Spin gap of the spin-1 tube as 
a function of the anisotropy parameter. Jm = 0.1J± and 
both open and periodic boundary conditions (OBC/PBC) are 
considered. 
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Figure 14: (color online) Spin gap of the spin-1 tube as a 
function of the anisotropy parameter (open symbols) or of 
the effective 2-leg spin-1/2 ladder (filled symbols). J« = 0.1J± 
and open boundary conditions are used. The infinite-size ex- 
trapolation was performed with the two largest ladder clus- 
ters. 



two biggest ladders (L = 64 and L = 128) which confirm 
that the spin gap has a large drop around a = 0.57. Ho- 
wever, at the critical point, the spin gap does not vanish 
and our data suggest an extremely small but finite value. 
This result would indicate a first-order phase transition, 
in agre ement with other numerical studies on the ladder 
systems^™ 



3. Von Neumann entropy 

In this part, we fix Ju/J± = 0.1 for which we have 
found two quantum phase transitions for a = 0.57 and 
a = 1.92, both with an extremely small but finite spin 
gap at the transition. Since these critical points seem to 
be very close to tricritical points (the correlation lengths 
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Figure 15: (color online) Plot of the Von Neumann entropy 
as a function of d(£) for different values of the anisotropy 
parameter and fixed J« = 0.1J±. 



at the critical points, which is valid on rather large length 
scales. Therefore, we believe that this spin tube can be 
tuned very close to a tricritical point separating a first- 
order transition line from a possible continuous transition 
line nearby. The gapless transition at the critical points 
could then be correlated with the presence of a non tri- 
vial topological term in the NLcrM. Note that there are 
no such critical points in the phase diagram of the two- leg 
ladder with S = 1 | 1Q | 49 1 , and that there is no topological 
term in the corresponding NLcrM as welP^. Of course, 
the link between the NLcrM and the critical point of the 
DMRG data still needs to be clarified with further inves- 
tigations, both analytically and numerically. 



VII. 



THE S = 2 CASE 



are very large at the transitions), another quantity of 
interest here is the von Neumann entropy of a finite seg- 
ment of the chain. It is defined by : 



S v n(£) = -Tr(pi\npt) 



(48) 



where pi = Trip is the reduced density matrix associated 
to a block of £ spins. This quantity is known to behave 
fundamentally differently for critical and non-critical sys- 
tem d 45 * 4 ^ . It saturates at a finite value when the system is 
non critical while it increases logarithmically for critical 
systems. The analytic expression of Si is given by : 



S v n(£) = 7 In 
6 



21 



sin 



9 



(49) 



where c is the central charge, and g is a constant. The 
Von Neumann entropy is represented on Figure [15] as 
a function of the conformal distance d(£) = s'm(ir£/L). 
For a — 0.5 and a = 0.6, the entropy converges to a 
finite value (the same results holds for all other values 
far enough from both critical points). On the contrary, 
for a = 0.57 and a = 1.92, it does not saturate but grows 
with the system size, which should be in favor of a gapless 
character of the two critical points. The entropy displays 
also a large periodic oscillation. Such an oscillation has 
already been observed in other critical spin chains 4 ^ and 
may be related to the existence of soft modes at k = and 
fc = 7r in the problenPSI. Thus, the decay in correlation 
function would not be simply algebraic at the critical 
point but the decaying function would be multiplied by 
an oscillatory factor. Because of the oscillations, it is hard 
to distinguish what is the best fit between c = § and 
c = 2. 

However, from our DMRG data on the gap, as well 
as the mapping to an effective ladder, we conclude to a 
very weak first-order transition, especially in the strong 
coupling limit Jy <C Jj_. This last scenario is supported 
by the bosonisation studies and the DMRG simulation 
of the effective Hamiltonian (44). Still, von Neumann en- 
tropy exhibits a critical behaviour with c = 3/2 or c = 2 



In this section, we consider the spin-2 case for which 
four quantum phase transitions are expected as a func- 
tion of the frustration a for a fixed J\\/J± <C 1. 



A. Effective model for the S — 2 tube 

In principle, one can apply simple perturbation theory 
in order to obtain an effective model valid close enough 
to any of the critical points for small J\\/J±- Here, we 
have two types of critical points : (i) close to a = 1/2 
or a = 3, and as recalled in Sec. |Hj a single triangle will 
have a quintet and a triplet as low-energy states ; (ii) 
close to a = 2/3 or a = 3/2, low-energy states consist in 
one triplet and one singlet. 

Although one can derive both kinds of effective mo- 
dels, case (i) does not allow to make analytical progress. 
On the contrary, for the second case, the effective model 
turns out to be formally the same as for the spin-1 tube 



(see part VIA), i.e. first-order degenerate perturbation 



results can be mapped onto a SU(2) spin-1/2 ladder that 
only contains 2-spin exchange interactions of the form 
given in Eq. (44). The effective exchanges are given as : 



2/3 : J ± = 3ev - 2, J\\ = 
7 

Jd= gJ||. 

3/2 : J ± = 3 - 2a, Jn = 



23 



-Jn 



(50) 



151 7, 7 39 7, 

^ J H> Jd -40 J " 



As explained in details in part VIA such a mapping ex- 
plains the occurence of a quantum phase transition when 
J± ~ 2Jrf, as well as giving insight on the order of the 
transition. 

Both quantum critical lines are plotted in Fig.j5] Mo- 
reover, for both cases we are in a regime where J&j ' J x ~ 
0.25 — 0.30, for whic h th ere is no consensus yet on the or- 
der of the transitiorPHS. Still, if the transitions are first- 
order, the gap at the transition should be very small, 
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Figure 16: (color online) Various spin gaps A a b for the spin-2 
tube as a function of a for J\\/J± = 0.1. Data are extrapolated 
to the thermodynamic limit from simulations with systems of 
length L = 8, 16, and 32. 
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Figure 17: (color online) String order for the spin-2 tube as 
a function of a for various lengths L with (a) J\\/J± =0.1 or 

(b) J||/J X = 1. 



which means that for most practical purpose, the sys- 
tem will appear critical on length scales smaller than the 
(large) correlation length. 



B. DMRG results for the spin 5 = 2 tube 

Simulations are done mostly with open boundary 
conditions (OBC) with system sizes up to 3 x 64, but 
also with periodic boundary conditions (PBC) on some 
cases. Typically, we keep up to 1600 states, which is suf- 
ficient to have a discarded weight smaller than 10~ 6 . 

We have computed several spectral gaps with various 
boundary conditions and various total spin sector, in or- 
der to avoid edge effects, but finite-size effects are rather 
large and spin gap values quite small so that no definite 
answer on the phase diagram can be obtained this way. 

A clearer signal is given by the string order parame- 
ter (see Eq. 47 1 which can distinguish between odd- J 



and even- J Haldane phases. Unfortunately, this quantity 
alone cannot distinguish between J = and J = 2 phase 
but we can rely on the presence/absence of edge states to 
distinguish these phases. In Fig. [TBJ we plot various spin 
gaps as a function of the frustration a. Obviously, if edge 
states carry a spin Sb, then the bulk spin gap is obtained 
between lowest levels of total spin (2Sb) and (25& + 1), 
A25 bi 2S b +i, i.e. once edge states have been polarized. In 
the thermodynamic limit, the two Sb edge states form 
(2Sb + l) 2 degenerate states. Our data are obtained on 
finite lengths up to L = 32 and we perform a 1/L linear 
extrapolation to get an estimate of the thermodynamic 
limit spin gaps. For a small J\\/J± = 0.1, our data shown 
in Fig. 16 indicate successively regions with Sb = 1, 1/2, 



and then in reverse order for increasing a, as had been 



conjectured initially in section III B 



In order to get more insight in this phase diagram, we 
have computed the string order parameter (see Eq. 47) 
for several exchange couplings. Data are shown in Fig. |17| 
and confirm the existence of even/odd Haldane S phases : 
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Figure 18: (color online) Numerical phase diagram for the 
spin-2 tube obtained from the spin gap and the string order 
parameter. 



a spin phase has a vanishing order parameter already 
on finite systems (see a close to 1) ; a spin 1 phase has 
a finite and positive string order parameter ; a spin 2 
phase should have a zero string order parameter but the 
finite-size scaling is rather slow, as is already known for 
the spin-2 Haldane chain for instance^. Note that the 
string order data are perfectly consistent with the edge 
states picture drawn from the spin gap data. This way, 
one can draw a quantitative phase diagram for the spin 
2 tube in Fig. |18[ , which confirms the qualitative picture 
conjectured in Fig. [5] 



VIII. 



DISCUSSION AND CONCLUSION 



In this paper we have shown that a simple model such 
as a three-leg quantum spin ladder can give rise to a very 
rich phase diagram. As it is now ubiquitous in quantum 
spin chains, integer and half integer spin must be treated 
separately. For half-integer spins, the Berry phase ana- 
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lysis of Section V C 1 points towards a quantum phase 



transition between a gapped spectrum and a degenerate 
ground state for a close to 1 and a gapless regime on each 
side of this phase, as it has indeed being observed for the 
5 = 1/2 case^. The semiclassical picture of this scenario 
is a phase transition separating a gapped isotropic co- 
planar phase and a pseudo-collinear gapless regime. Not 
surprisingly the difference in behavior at large scales is 
dictated by the Berry phase terms present in the action. 

For integer spins the situation is even more interesting. 
If we consider the strong coupling regime J± ^> Jii it 
makes no doubt that 25 quantum phase transitions are 
expected for a spin 5 tube when varying the anisotropy 
parameter a. These phase transitions separate gapped 
phases, and this scenario is reminiscent of what happens 
in dimerized spin chains when varying the dimerization 
parameter (see for example I. Affleck's lectures^) . 

These phase transitions can be understood in terms of 
spontaneous breaking of the hidden Z 2 x Z 2 symmetry. 
The broken hidden Z 2 x Z 2 symmetry can be detected 
by a string order parameter or by edge states with half- 
integer spin. While a simple generalization of the string 
order parameter to the tube vanishes in the weak rung 
coupling limit, an alternative string order parameter re- 
mains finite in the same limit for odd spin 5. Additio- 
nally, in some regions of the phase diagram, there exist 
"phases" with integer-spin edge states. Although they ap- 
pear to be non-trivial phases, they can be adiabatically 
connected to a trivial phase with no edge state. This is 
consistent with an unbroken hidden Z 2 x Z 2 symmetry 
in these "phases" . 

More insight into the phase diagram can be obtained 
with the recent discussi on conc erning the characteriza- 
tion of the Haldane phase^J22152j rpj^ jjaldane phase can 
remain a distinct phase separated from a trivial phase 
by a quantum phase transition, even when the hidden 
Z 2 x Z 2 symmetry is not well-defined and the string or- 
der parameter is not useful. It turns out that the Haldane 
phase is a topological phase protected by any one of the 
following three symmetries : 1) global D 2 (= Z 2 x Z 2 ) 
symmetry of 7r-rotation about x,y and z axes, 2) time- 
reversal symmetry (for Sj — > —Sj), and 3) lattice inver- 
sion symmetry about a bond center (link parity). The 
hidden Z 2 x Z 2 symmetry is well-defined only with the 
symmetry 1) above. Most generally, the Haldane phase 
is characterized by an exact double degeneracy of the 
entanglement spectrum. 

In this paper, we limited our discussion to the tube 
with SU(2) symmetry of spin rotation and all the sym- 
metries l)-3) listed above. Within this limitation, the 
hidden Z 2 x Z 2 symmetry can be used to characterize 
the nontrivial phases, which have edge states with a half- 
integer spin Si,. On the other hand, we expect that the 
phases with the broken hidden Z 2 x Z 2 symmetry cor- 
respond to the generalized Haldane phase with an exact 
double degeneracy of the entanglement spectrum. It is 
protected by either of the symmetries 2) or 3), even 
when the symmetry 1) is explicitly broken and the hid- 



den Z 2 x Z 2 symmetry is ill-defined. In particular, as long 
as the lattice inversion symmetry abount a bond center 
is preserved, the generalized Haldane phase is protected 
as a distinct topological phase. This protection may be 
roughly understood by the intrinsic odd parity with res- 
pect to the lattice inversion associated to an odd number 
of valence bonds between the neighboring runga^. 

The above general analysis implies the existence of a 
quantum phase transition between a "topological" phase 
and a "trivial" phase. However, it does not tell the order 
of the transition or its universality class. 

For dimerized chains, the NLerM approach shows that 
the critical points correspond to an effective half-integer 
chain which is described by an SU{2)\ WZNW model. 
Our analysis of the NLerM in the triangular spin tube 
has shown that the situation is different here indica- 
ting that we must expect phase transition of a different 
kind. Arises then the question of whether these transi- 
tion points are expected to be first order or more interes- 
ting critical theories as for example higher levels SU(2) 
WZNW models. 

The case of the 5 = 1 ladder has proven to be a very 
interesting and instructive example. The low energy be- 
havior of this system can be shown to be equivalent to 
a two-leg spin 1/2 ladder. This ladder system has two 
obvious extreme regimes corresponding to a collection 
of singlet states (strong antiferromagnetic couplings bet- 
ween the chains) and a Haldane phase of an effective 
spin 1 chain (strong ferromagnetic coupling between the 
chains). The most recent bosonization analysis^ indicates 
that the transition between these two regimes can be ei- 
ther first order, or a couple of (gapless) lines surroun- 
ding a dimerized phase. These gapless lines have central 
charges c = 1/2 and c = 3/2, this last one corresponding 
to a 5[/(2) 2 WZNW model. We have performed exten- 
sive DMRG computations on this system. The analysis of 
the spectral gap and the von Neumann entropy tend to 
indicate a weak first order transition for our system, but 
in any case the close proximity to the tricritical point. 
This allows us to speculate that by introducing further 
microscopic parameters, as for example second neighbors 
interactions, the transition can be made second order, 
but this issue is beyond the scope of the present work. 
This result is also encouraging for analyzing the nature of 
the transition for higher spins with both numerical and 
novel analytical techniques. 

One important point is that many of the results obtai- 
ned here generalize to ladders with a higher odd number 
of legs displayed with periodic boundary conditions (so 
in a frustrating manner) . Of course frustration becomes 
weaker as one increase the number of legs. In this sense 
the 3-leg ladder is a representative of a family of quasi 
one-dimensional systems where frustration plays a crucial 
role in the emergence of an interesting physics. 
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